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Abstract 

We analyze statistical consequences of a conjecture that there ex- 
ists a fundamental (indivisible) quant of time. We study particle dy- 
namics with discrete time. We show that a quantum-like interference 
pattern could appear as a statistical effect for deterministic particles, 
i.e. particles that have trajectories and obey deterministic dynamical 
equations, if one introduces a discrete time. As a demonstration of 
this concept we consider particle scattering on a screen with a slit. 
We study how resulting interference picture depends on the parame- 
ters of the model. The resulting interference picture has a nontrivial 
minimum-maximum distribution which vanishes, as the time discrete- 
ness parameter goes to zero. This picture is qualitatively the same 
as one obtained in quantum experiments. The picture includes some 
interesting nonclassical properties such as a 'black' region behind the 
center of the slit. 
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2002-00279 (Network on Quantum Probability and Applications) and Profile Math. Mod- 
elling in Physics and Cogn. Sc. of Vaxjo University. 
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1 Introduction 



In papersPQ 13 "we presented the conjecture: 

there exists a fundamental quant of time r. 

In fact we came to this conjecture by finding (by pure occasion) that sta- 
tistical histograms obtained for classical systems in the one or two slit ex- 
periments can have the form of quantum-like interference observed in ex- 
periments with photons or electrons. We started to speculate that essential 
features of quantum statistics (at least some of them) could be obtained as 
a consequence of existence of time quant. 

Basically speaking the idea of discretization is quite natural for physics. 
Indeed discretization of energy was exploited by Plank-Einstein ideas in order 
to explain such quantum phenomena as black body radiation, photo-electric 
effect, etc. On the other hand from the point of view of energy-time un- 
certainty relations discretization of time gives certain constraints to a (mea- 
sured) energy of the system. 

On the other hand by introducing discretized time one effectively changes 
the theoretical structure of the space-time on small distances^. From this 
point of view it is rather natural to assume that a value of the discreteness 
parameter is of the order of Plank time tpi ^ 5 - 10~'*^(sec). One could think 
of discretized time as a lack of information of where the particle is during 
the discretization period. 

In this paper we study deterministic model for particle scattering on 
a screen with a single slit. Resulting interference picture has a nontriv- 
ial minimum-maximum distribution which is qualitatively the same as one 
obtained in quantum experiments. The interference pattern appears as a sta- 
tistical effect of particles which have trajectories^. The interference picture 
has no relation to 'self-interference' of particles, no wave-structure is involved 
into considerations. The basic source of interference is the discrete time scale 
used in our mathematical model: instead of Newton's differential equations 
with continuous time evolution, we consider difference equations with dis- 
crete time evolution. Interference effect disappears as the time discreteness 
parameter goes to zero. The picture includes some interesting nonclassical 
properties such as a 'black' region behind the center of the slit^. 

^In particular, this idea was exploited by p-adic space-time alternative [HI 1101 1111 112j 
^Cf. with Bohmian mechanics |lf)l ITT) 

•^We would like to thank G. Emch for an interesting question he asked on the 'Interna- 
tional Conference: Reconsideration of Foundations-2' held in Vaxjo - whether it is possible 
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The paper is organized as follows. In the next section we provide basic 
ideas of the discrete time formalism and compare it with construction of 
quantum mechanics. In section El we write dynamical equations with discrete 
time and in section |3] we present them in intuitive form. In section El we 
apply the formalism to a particle scattering on a screen with a slit and study 
appearance and properties of the resulting interference picture. Finally, in 
the appendix we provide a detailed description of the numerical simulation 
performed, we list a Mathematica program which mirrors the original parallel 
C++ programm which we used to perform actual simulation. 

We remind that recently it was suggested jH El HI IB] that quantum 
mechanical interference rules could be explained using contextual transfor- 
mations. This approach makes one think of a construction of physically 
reasonable model which produces interference resting on the classical rule 
of addition of probabilistic alternatives (extended with a notion of context) . 
Our quantum-time approach could be considered as a candidate for such a 
model. 

2 Basic Steps of the Approach 

According to the principles of quantum mechanics in order to compute mea- 
surable quantities, in the simplest case, one has to perform the following 
steps |13| in]. First, the form of classical and quantum equations of motion re- 
mains the same (we are talking about Heisinberg's representation), although 
the dynamical variables are now noncommutative. Next, one has to add to 
the original initial conditions an extra one stating commutative properties of 
the canonical variables, like [x(0),p(0)] = ih. At the last step one says that 
measurable quantities are eigenvalues of the dynamical variables. Here the 
second step introduces the Plank constant which distinguishes classical and 
quantum worlds. The noncommutative property of dynamical variables is the 
root of the fact that quantum particles can not have trajectories and intro- 
duction of quantum state implies the nondeterministic behavior of quantum 
systems. 

In this note we tried to achieve the same results for measurable quantities 
as in quantum mechanics (at least in the part vastly supported by experi- 
ments), but trying to keep trajectories and deterministic behavior of parti- 

to reproduce Fresnel-like phenomenon of a black region behind the center of the sUt in the 
discrete time formahsm - our answer now is yes (see section |S1 for details). 
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cles. We follow steps similar to the above ones, but we keep the measurable 
quantities being real functions. Nevertheless, we have had to introduce the 
time discreteness parameter r. We start from classical equations of motion 
with discrete time and postulate the canonical variables to be measurable. 



3 Discrete Time Dynamics 

In both classical and quantum mechanics a dynamical function A = A[p, q) 
evolves according to the following well known equation 

DtA=[H,A] (1) 

where H = H{p, q) is a Hamiltonian of the system and in the right hand side 
is a Poisson bracket, either classical or quantum^. 

In both classical and quantum dynamics the left hand side of (Q) is the 
same, it contains a continuous time derivative 

DtA = — 
dt 

As it was mentioned in the previous section we are interested in con- 
struction dynamics with discrete time. This is done with the help of discrete 
derivative which is postulated to be 

Dl^^A = -[Ait + T)-A{t)l 

T 

where r is the discreteness parameter. This parameter is finite and is treated 
in the same way as Plank constant in quantum mechanical formalism. In 
particular if r is small relative to dimentions of the system then classical 
approximation with continuous derivative should work well (although this 
could not be the case all the time in the same sense as there are examples 
when quantum formalism is reasonable even for macroscopic systems, for 
example in superfluidity). 

^Poisson bracket in classical mechanics could be presented as 

[A, B] (2) 
and in quantum mechanics it is a commutator [A, B] — j^r^AB — BA) 
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Thus the discrete time dynamical equation is postulated to be 



where A{p, q) is a real- valued function of real- valued dynamical variables and 
in the right hand side there is the classical Poisson bracket (see 

Please note that in our model the coordinate space is left continuous. 



4 Discrete Time in Newton's Equations 

In this section we provide an intuitive description of the discretization pro- 
cedure described above. Consider a well known Newton's equation 



mr 



(3) 



This equation gives the same dynamics as (H}-©- Let us now modify this 
equation, introducing the time discreteness parameter r. We rewrite the 
second order differential equation (jH)) as a system of first order differential 
equations, we have 

F = mi) , , 

4 

V = r 

In the system (0]) the derivatives assume continuousness of time. Let us now 
introduce the discreteness parameter r. We have 

v{t + r)-v{t) 
t = m 

r(t + r) - r{t) ^ ' 



T 



In the limit of r ^ © is equivalent to Q and (jD). 

Now let us rewrite the system (0) in a directly computable way 



v(t + t) = v(t) + Fr/m 

r{t + r) = r{t) + v{t)T ^ ^ 



where F = F{r{t),v{t),t). 
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Figure 1: Particle scattering on a screen with a slit. Particles are emitted 
from the source e pass through a slit in screen Si and gather on screen S2. 
There is nonzero force field in the region S near the slit. 



5 Interference Pattern in a Single Slit Scat- 
tering 

Consider the following experimental setup (see Fig^. A particle source e is 
located in front of the center of a slit in a screen Si. Near the slit there is a 
nonzero force field which affects particles, 

Fix,y) ^ 0, for {x,y) E S 

{x and y are coordinates along horizontal and vertical axes respectively, the 
axes origin is in the center of the slit). Particles pass through the slit and 
concentrate on a second screen 5*2. We study a particle density on the screen 

Let us start from the simplest case when the force field is constant and 
perpendicular to the screens, F{x,y) = Fq, for x ^ and otherwise. 
Let the source e be point-like emitting particles with constant velocity vq 
under random evenly distributed angles a G [0, 2n). In this case trajectories 
of particles emitted by e (in discrete time dynamics) in the region a; < (i.e. 
before the screen 5*1) form concentric circles originating from e with the radii 
f^n = VoTTi, n = 1, 2, . . . (Fig|21). We get the circles - which are trajectories 
of many particles emitted under close angles - in the region where there is 
no force field and particles move along strait lines, exactly as in classical 
dynamics. Let be the points where these circles enter the region 6. The 
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Figure 2: Particles emitted from the source e enter the region 5. Points a±i, 
a±2, a±3, . . . are origins of deviation which forms interference pattern. 



distance between the center of the sht and is given by 



where d is a distance between e and Si, and no is the largest integral value 
not greater than a fraction d/voT, 



This rather simple setup already produces interesting nontrivial interference 
pattern (FiglSl). 

The points are the origins of deviation from 'classical' (i.e. continuous 
time) trajectories. This deviation forms the interference pattern. To argu- 
ment this let ipi be an angle between the horizontal axis and a line connecting 
e and Oj. Particles emitted under angles less than ipi {i^ 1) become affected 
by the force field in the region 5 one step earlier than those emitted under 
angles grater than ^pi. As a result there appears a 'fork' - even very close 
trajectories but with the angles above and below ipi become separated. One 
could get the points of minima of the interference pattern by following the 
trajectory along the line e ^ and the parabolic curve (movement in a 
constant field) in the region x > until the trajectory hits the screen 5*2. 

The case described above is not completely physically reasonable, indeed 
is requires a force field everywhere in the half-plane x ^ 0, thus it could be 
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Figure 3: Interference pattern (brighter area corresponds to the higher par- 
ticle density) and particle trajectories in a single-slit scattering computed in 
a discrete time formalism. 

considered only as a simplified model which still demonstrates some inter- 
esting properties of discrete time dynamics. An example of such a property 
is a Fresnel-like phenomenon of a black region behind the center of the slit 
which is discussed below. 

One could see from dZI)-® that if 

d = vork, , where k = 1,2, . . . 

then a+i = a_i, and thus the particle trajectories emitted under positive and 
negative become separated. This forms a black region in the center of the 
screen 5*2 just behind the slit (FigE)). 

Although we succeeded with a model of a field localized near the slit. We 
took Gaussian rapidly changing force field of the form 

F{x,y) = e,Foe-''^\ 

It is very stimulating that such setup also produces nontrivial pattern, since 
it means that one could start thinking of a physical nature of this force field. 

We have also considered the case of non-point-like Gaussian source dis- 
tributed along the vertical line. This produces smoother picture, which could 
be obtained by convolving the pattern on Fig|21with a Gaussian kernel. 
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Figure 4: Fresnel-like phenomenon of a black region behind the center of the 
sht in the discrete time formahsm. 
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Numerical Simulation 



Here we provide a Mathematica program for computation of the histogram. 
Please note that although this program is approximately 10^ times slower 
than a more powerful parallel C++ program which we used to actually per- 
form computations it reproduces the same results (in the simplest case where 
the computation time allowed us to wait the result) and could be seen as a 
detailed description of the numerical experiment which was performed as a 
demonstration of the concept. 
ComputeTraj ectory ^Compile [{{a, _Real} } , 
With[{ 

d = 5,1 = 10,6 = 0.5, R = 5, {* Frame parameters *) 
q = -1, (* Charge *) 
T = 0.025, (* Discreteness parameter *) 
= 12 (* Initial velocity *) 

}, 

Module [{ 

X — —1.0 * d,y = 0.0, Vx = V * Cos[a], Vy = v * Sin[a], Fxv = 0.0, Fyv = 
0.0}, 

While[x < I, 

Fxv = If[^&s[x] ^ 5, 27rq, 0]; 

Fyv = 0; 
x+ = Vx * r; 
y+^Vy* t; 
Vx+ — Fxv * r; 
Vy+ — Fyv * r; 

]; 

Module[{ 

xPrev — X — {Vx — Fxv * r)r, 
yPrev — y — {Vy — Fyv * t)t}, 

(y-yPrev)^E^^+yPrev]]]]; 
ComputeHistogramm[r _ : 0.1] := Module[ 

{b, i = I, iPrev, t = {}, h}, 

For[b = tbl[[l]] +r,b< tbl[[-l]], b+ = r, 
iPrev — i; 

While[tbl[\i]] <b,i + +]; 
{t,h[{b,i-iPrev}]} 
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Flatten[t]/ .h[a_] — > a 

]; 

Timing[tbl = Table[ComputeTrajectory[a°], a, —49, 49, 0.01]; ] 
hst = ComputeHistogramm[0.025]; 
ListPlot[hst, Plot Joined — > True, PlotRange — > All]; 
ListContourPlot[Transpose[hint / .X _, P_ — > P, P], 

ContourLines — > False, Contours — > 100, FrameTicks — > False] 
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